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Abstract. An MHD theory is developed for the stand- 
off distance a s of the bow shock and the thickness A ma 
of the magnetosheath, using the empirical Spreiter et 
al. relation A mj = kX and the MHD density ratio X 
across the shock. The theory includes as special cases 
the well-known gasdynamic theory and associated phe- 
nomenological MHD-like models for A ma and a s . In 
general, however, MHD effects produce major differ- 
ences from previous models, especially at low Alfven 
(Ma) and sonic (Ms) Mach numbers. The magnetic 
field orientation, M^, Ms, and the ratio of specific heats 
7 are all important variables of the theory. In contrast, 
the fast mode Mach number need play no direct role. 
Three principal conclusions are reached. First, the gas- 
dynamic and phenomenological models miss important 
dependances on field orientation and Ms and generally 
provide poor approximations to the MHD results. Sec- 
ond, changes in field orientation and Ms are predicted 
to cause factor of ~ 4 changes in A ma at low Ma- These 
effects should be important when predicting the shock’s 
location or calculating 7 from observations. Third, us- 
ing Spreiter et aids value for k in the MHD theory leads 
to maximum a s values at low Ma and nominal Ms that 
are much smaller than observations and MHD simula- 
tions require. Resolving this problem requires either 
the modified Spreiter-like relation and larger k found 
in recent MHD simulations and/or a breakdown in the 
Spreiter-like relation at very low Ma- 

1. Introduction 

The location of Earth’s bow shock has been actively 
researched since its prediction and discovery. Subjects 
of particular interest include the shock’s farthest ex- 
tent sunwards (known as the standoff distance a 9 ) and 
the thickness A m , of the magnetosheath region separat- 
ing the shock from the magnetopause, due to their im- 
portance in understanding foreshock observations, so- 
lar wind-magnetosphere interactions, and the ratio of 
specific heats 7 for the plasma. Figure 1 defines a s , 
A ma , and the magnetopause standoff distance a mp in 
the X-Y-Z coordinate system formed by rotating the 
GSE system so that the shock is symmetric about the 
solar wind’s velocity vector relative to Earth. Clearly 
a 5 — a mp + A ms . Balancing the solar wind ram pressure 
P — p sw v 2 sw and the magnetostatic pressure of Earth’s 
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dipole magnetic field requires a mp = K P _1 ^ 6 , where A 
is a slowly varying function of the IMF B z component, 
the ring current, the magnetopause current system and 
drag effects for the solar wind-magnetosphere system 
[Formisano et al , 1971; Slavin and Holzer, 1978; Far- 
ris et al., 1991; Sibeck et al, 1991]. Spreiter et al. [1966, 
and references therein] developed the first detailed the- 
oretical model for a s and A m s > 


= I<P- 1/6 [1 + 1.1 


(7- l)M 2 + 2 


(1) 


v (7 + 1)A/ 2 

The ratio A m ,/a mp is given by the entire second term, 
with the number 1.1 depending on the obstacle’s shape. 
This equation follows from the empirical linear relation 


Ait = kX = k P — (2) 

&mp Pd 

obtained from gasdynamic simulations for M — Ms ~ 5 
[Spreiter et al, 1966] with the ratio p sw /pd specified 
subsequently by the jump conditions for a gasdynamic 
shock. (Here pd is the mass density downstream from 
the shock.) An analytic explanation for (2) remains 
unavailable. Spreiter et al generally identified the 
(sonic) Mach number M with the Alfven Mach num- 
ber Ma — v sw /va for arbitrary magnetic field orien- 
tations and M\ » Mj » 1 (a pseudo Mach number 
was defined instead for aligned flows with v sw || B au; ). 
This theory is therefore intrinsically gasdynamic with 
the subsequent phenomenological replacement of the 
sonic Mach number Ms = v sw /c 3 by Ma- Spreiter 
et al emphasized the theory’s expected limitations at 
low Mach numbers. More recently Russell [1985] sug- 
gested that the proper replacement for M in (1) is the 
fast magnetosonic Mach number M ms , since the bow 
shock is a fast mode shock. Again, however, this is a 
phenomenological replacement in a gasdynamic result. 

Earth’s bow shock is indeed a fast mode shock, not a 
gasdynamic shock, and so MHD theory is a priori more 
appropriate. Spreiter et al y s gasdynamic equation and 
its phenomenological variants therefore need to be re- 
considered and a MHD version of (1) derived more rig- 
orously. Other motivations for studying the bow shock’s 
location include the finding that (1) with M replaced 
by Ma or M ms predicts a s values that are too small for 
Ma & Mms ~ 1 — 3 [Russell and Zhang, 1992; Cairns et 
al, 1994] and the scattered values for 7 extracted from 
the measured A ma via (1) [Fairfield, 1971; Zhuang and 
Russell, 1981; Farris et al, 1991]. Lastly, MHD predic- 
tions for a s and A m , are needed now that global MHD 
simulation codes are available [e.g., Cairns and Lyon, 
1994] for studying the Earth-solar wind interaction. 
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Figure 1 . Definitions of a s , a mp , A ms , the angle 0, 
and the X-Y-Z coordinate system. 

This paper addresses the basis of previous gasdy- 
namic and phenomenological MHD-like models for a 3 
and A ma , the importance of MHD effects, and current 
attempts to explain unusually distant shock locations 
and the plasma’s ratio of specific heats 7. The simplest 
MHD theory for a 3 and A ma is constructed: Spreiter et 
ai ’s empirical relation between A ma and X is retained 
and MHD theory is used to specify the density ratio 
X across the shock. This MHD theory (Sections 2 and 
3) includes the gasdynamic and related phenomenolog- 
ical theories as special cases but in general shows differ* 
ent theoretical dependances. In fact, the MHD theory 
predicts that A ma and a 3 should depend strongly on 
the magnetic field orientation, Ma and Ms (and 7). 
Zhuang and Russell [1981] previously developed an ap- 
proximate but not self-consistent MHD expression for 
X and an unrelated calculation for A ma when Ma & 
Ms 1. This paper’s theory extends Zhuang and Rus- 
sell's work on X , by retaining all contributions to X 
and not assuming high Ma & Ms flows, and merges it 
with Spreiter et ai s empirical approach. Quantitative 
comparisons are made between the various theories in 
Section 3. A discussion and the conclusions are pre- 
sented in Sections 4 and 5, respectively. 


2. Analytic theory 


Spreiter et ai [1966] used gasdynamic simulations to 
show that the linear relation (2) between A ma and X 
holds for Ms ~ 5. The jump conditions for a gasdy- 
namic shock lead to a quadratic equation for X whence 


(7 - 1)M% + 2 

(7 + 1)MJ 


(3) 


thereby leading to (1) via (2). The transition to the 
magnetised solar wind was then attempted by phe- 
nomenologically replacing Ms by Ma [Spreiter et ai, 
1966] or M m3 [ Russell, 1985]. The limiting values 
of A ma and a 3 for 7 = 5/3 are then: A ma /a mp — > 
Ll/4 and a 3 /a mp — ► 1.275 as Ma and M ms — > 00; 
A ms /amp -4 1.1 and a s /a mp 2.1 as M A and M ms — > 
1. Note that the gasdynamic and phenomenological re- 
sults have no explicit dependences on magnetic field ori- 
entation and Ms except through M ma . 

The theory developed here assumes that (2) remains 
valid, as supported by Cairns and Lyon's [1994] MHD 
simulations for Ma ~ 1.5, and uses MHD theory to 
specify the density jump X. Factoring out the solu- 


tion X — 1, the MHD jump conditions lead to a cubic 
equation for X [e.g., Zhuang and Russell , 1981]: 

AX 3 + BX 2 + CX + D = 0 (4) 

A = (7 + 1 )M\ 

-B = (7 - + (7 + 2) cos 2 9M\ + (7 + 7 0)M\ 

C = (7 - 2 + 7 cos 2 9)M\ + (7 + 1 + 27/?) cos 2 0M\ 
— D — (7 — 1) cos 2 6M\ + 7/? cos 4 6 . 

The new variables introduced are 6 C [0,90°], the angle 
between v au , and B au , (the shock normal is antiparal- 
lel to v au >), and the upstream plasma /? is defined by 
7 f} = 2c 2 $ jv 2 A = 2M\jM\. The natural Mach numbers 
in the theory are then Ma and M5; M ms need play 
no role. In comparison, neither 6 nor /? play roles in 
the gasdynamic or phenomenological theories (except 
through Mms)- The standard cubic analysis provides 
general solutions to (4), although these typically pro- 
vide little insight. Consider, however, the special cases 
of parallel (0 = 0°) and perpendicular (6 = 90°) flows. 

The cubic is easily factored when B au , is parallel to 
the shock normal (cos0 = 1). Ignoring the two ‘switch- 
on’ shock solutions X — Af^ 2 , (2) and (4) yield 

_ , . , (7-1)Mj+7/? _ . , (7~l)Mj+2 

amp" (7 + l)M 2 + (7+1 )Mj ' 

(5) 

The rightmost form reveals that the gasdynamic ex- 
pressions for X and a s are recovered, cf. (1) and (3), 
as expected since the magnetic field essentially drops 
out of the problem in the parallel case. This equa- 
tion implies two important results for 0 = 0°. First, 
phenomenological replacement of Ms by Ma or M rns 
(= M a for 0 — 0°) in (1), as proposed by Spreiter et 
ai and Russell, is incorrect except in the special case 
7 /? = 2. The phenomenological theories are therefore 
restricted special cases of the MHD theory. Second, 
A ma , a 3 and A" are independent of Ma and M ms , and 
depend only on 7 and Ms (with the intuitive caveats 
that Ma & Ms > 1). This is a major difference from (1) 
with the replacements M — > Ma or M ms , which predict 
A rn 9 /a mp k as Ma — > 1 instead of the correct result 
A rns/amp “> k / 4 for 7 = 5/3 and M s » 1. 

For perpendicular flows (cos# = 0) the cubic equation 
for X collapses to a quadratic, whence 

( 6 ) 

with A = (7 — 1) + 7 /M\ + 2/ Mg. For 7 < 2 only 
the solution a 3 = a 3 + is relevant (only X = X+ > 
0). In general, (6) is not equivalent to the gasdynamic 
solution (1), the Spreiter et ai and Russell variants of 
(1) with M — > Ma and M ms , respectively, or (5)’s MHD 
solution for 0 — 0°. However, in the special case 7 = 2 
(6) reduces to Russell's form for arbitrary f3 (writing 
Ma and Ms in terms of M ms and j3) and to Spreiter 
et ai ’s form in the limit fi — > oc. Furthermore, the 
gasdynamic result (1) is recovered in the limit Ma —► 
oo (or f3 — ► 0) for arbitrary 7. (This is true for all 
9 .) Equation (6) also implies three corollary results. 



First, a,, A m3 and X vary strongly with the angle 0. 
Second, in general the solution depends intrinsically on 
both Ma and A/5; while (6) can easily be rewritten in 
terms of A/ m3 and /?, only for 7 = 2 does the explicit 
1 3 dependance disappear and the solution depend solely 
on Mms . Last, the behavior as Ma — ^ 1 for 9 — 90° 
differs greatly from the MHD result (5) for 0 = 0: for 
7 = 5/3, (5) implies X — > 1/4 for 0 = 0 (Ms = oo) 
while (6) implies X+ — > 1 for ft = 0 and X+ —>1.7 for 
7/? — 1. Thus, factor of ~ 4 changes in A m< should 
exist at low Ma for different 0 and Ms . 

3. Numerical Analyses 

When the magnetic field is neither parallel nor per- 
pendicular to the flow direction and shock normal, (4)’s 
solutions are most instructive when presented graphi- 
cally. Figure 2 shows the MHD theory’s predictions for 
the ratio a s /a mp = 1 + 1.1X, where Spreiter et al.’s 
empirical value for k in (2) is used (see Section 4) and 
(4) is solved numerically, as a function of Ma for 0 = 0, 
45, and 90° and Ms = 8. Figure 3 is an analogous plot 
for Ms = 2. The strong dependences on 0, M a and Ms 
implied by (4) - (6) are clearly evident. Note that the 
0 = 0° curves are all flat, and so independent of Ma 
as in (5), with a level that depends on A/5. The curves 
for 0 = 45° lie below the 0 = 90° curves. Indeed, it 
may be shown analytically that the maximum allowed 
values of a s (and X) for given Ms occur at 0 = 90° and 
Ma — 1 + , with smaller Ms leading to larger X and a 3 . 

Figure 4 shows how X and a s depend on 7: for a 
given 9 and A/5, a larger 7 leads to a higher curve for 
a 3 versus M\- However, it can be shown analytically 
from (4) that the maximum values of X and a s are 
independent of 7: curves for different 7 therefore all 
converge to the same maximum at Ma — 1 and 0 = 90°. 

Quantitative comparisons between the MHD theory 
developed here, Spreiter et al. ’s model given by (1) with 
M = M a y and Russell's model are shown for various 0 
and M a in Figure 5. For all 0 the gasdynamic theory 
with M = Ms coincides with the 0 = 0° MHD solution. 
In general significant differences between the MHD pre- 
dictions and the phenomenological models are apparent, 
except in the limited region Ma ~ 1.5 and 0 ~ 30° (de- 
creasing Ms shrinks this region further). Note that 
the phenomenological models have no 0 or (direct) A/5 



Ma 


Figure 2. Ratio a s /a mp as a function of M A and 0 for 
7 — 5/3, Ms — 8 and k = 1.1. The ratio is independent 
of Ma for 0 = 0° and is maximum at 0 = 90° & Ma — 1. 



Figure 3. Similar to Figure 2 but with A/5 = 2. Note 
that a, increases as Ms decreases (for all 0). 


dependences, whereas the MHD theory shows these de- 
pendences to be very important. The maximum a, val- 
ues predicted by the MHD and phenomenological mod- 
els (M a ~ 1) are, however, almost identical for large 
Ms ~ 5. When A/5 is small the MHD theory predicts 
larger a $ values. However, the maximum difference in 
a 5 between the MHD and ‘gasdynamic’ theories remains 
less than 50% for A/5 > 1. 

4. Discussion 

Differences of 50% or more in a s and 400% in A m3 , 
due to 0 and Ms effects when Ma ~ 5, are easily dis- 
cerned in Figures 2-5. Differences occur both between 
the MHD, gasdynamic and phenomenological theories 
and between different parameter sets for the MHD the- 
ory. That 0 effects should be very important, as well as 
M a and Ms variations, in determining a s and A m5 is 
one of this paper’s important predictions. Refinements 
to the present ‘local’ MHD theory will undoubtedly oc- 
cur when global MHD effects are considered. 

The above results argue that values for 7 derived from 
measurements of A m3 [Fairfield, 1971; Zhuang and Rus- 
sell, 1981; Farris et al., 1991] should depend on the for- 
mulae used for A m3 and on whether Ma, 0 and A/5 
variations are all considered. Finite Ma, 9 and A/5 ef- 
fects may explain the wide scatter in the published 7 
values. Analyses to measure 7 should be redone using 
(4)’s explicit MHD solutions (or successors thereof). 

Near 1 AU the solar wind speed and electron temper- 
ature vary relatively little, whence A/5 lies within a fac- 
tor of 2 of the nominal value Ms ~ 7. For Ma ~ 1 — 3, 



Figure 4. Ratio a 3 /a m « versus Ma for 7 = 5/3 (full 
lines) and 7 = 2 (dashed), A/5 = 8, and 0 = 0 & 90°. 
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Figure 5. Comparison between the MHD theory given 
by (2) and (4) (full lines) for various 6, and the phe- 
nomenological Spreiter et al. (dashed) and Russell (dot- 
ted for 0 — 90°) models for Ms = 5 and 7 = 5/3. For 
arbitrary 9 the Russell curve lies between the dotted 
and dashed curves. 

then, Figures 2-5 predict that a s will be within 10% 
of (l)’s prediction (for M ~ Ma)- This paper’s MHD 
theory, using Spreiter et al.’s empirical value for k in 
(2) and MHD theory to determine A\ can therefore 
not explain unusually distant bow shock observations at 
low Ma [Russell and Zhang , 1992; Cairns et al. } 1994]. 
Cairns and Lyon’s [1994] MHD simulations show, how- 
ever, that intrinsically MHD effects modify (2) and in- 
crease k by a factor ~ 3 over Spreiter et al. \ s value. 
Incorporating the modifications into this paper’s MHD 
theory they can explain the large a s values at low Ma in 
the simulations. Comparisons with observational data 
still need to be done. It remains possible that the MHD 
variant of (2) breaks down at very low M a ■ 

5* Conclusions 

The analyses above show that more detailed consid- 
eration of MHD effects results in major differences from 
the well-known gasdynamic theory, and its phenomeno- 
logical variants, for A m5 and a s . The MHD model de- 
pends on Ma ^ Ms , 9, and 7, all of which have important 
effects; in comparison the gasdynamic and phenomeno- 
logical models involve only a single (differing) Mach 
number and no explicit 9 effects (save through M ms ). 
The MHD theory reduces to the gasdynamic theory in 
the special cases 9 — 0° or Ma Ms T The Spre- 
iter et al. [1966] and Russell [1985] phenomenological 
theories reappear as the special cases of the MHD the- 
ory for (i) 9 = 0 with 7/? = 2, and (ii) 0 = 90° with 
7 = 2 and /? = 0 or arbitrary /?, respectively. The MHD 
theory predicts that A m , and a s should depend strongly 
on 9 , Ma and Ms when Ma and/or Ms ~ 5, with the 
functional form of the Ma and Ms dependence varying 


with 0.1n particular, varying 9 from 0° to 90° at low M a 
should cause A mj to vary by ~ 400%. Since Ma and 
Ms are frequently of this order out to 1 AU, the pre- 
dicted MHD effects should have widespread applicabil- 
ity; for instance, in calculations of 7 from measurements 
of A mj . It is found that the MHD theory with Spre- 
iter et al. *s value for k cannot explain observations and 
MHD simulations of distant, very low Ma bow shocks. 
Cairns and Lyon [1994] can explain the simulation re- 
sults by inserting the larger, MHD value for k into the 
present paper’s MHD theory. It remains possible, nev- 
ertheless, that a nonlinear model for A mj = A ma (A) is 
necessary at very low M a - 
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